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Abstract 

Since present Monte Carlo algorithms for lattice QCD may become 
trapped in a fixed topological charge sector, it is important to understand 
the effect of calculating at fixed topology. In this work, we show that 
although the restriction to a fixed topological sector becomes irrelevant in 
the infinite volume limit, it gives rise to characteristic finite size effects 
due to contributions from all 6'-vacua. We calculate these effects and 
show how to extract physical results from numerical data obtained at 
fixed topology. 



1 



1 Introduction 



The numerical solution of lattice QCD is a notoriously hard problem. In particular, 
in order to obtain physical results, the numerical data must be extrapolated to 
the chiral, continuum, and infinite volume limits. Fortunately, the dependence of 
physical quantities on the quark masses, the lattice spacing, or the finite box size 
can often be understood analytically, and formulae can be derived that allow us to 
extrapolate reliably to the physical limit. Here we derive such formulae to correct 
for finite size effects that are due to fixed topology. 

The configuration space of QCD decomposes into topological sectors. With 
boundary conditions that are periodic up to gauge transformations the topologi- 
cal charge Q is an integer. In the continuum theory different topological sectors 
are separated by infinite action barriers. In a lattice theory the action barriers de- 
pend on the choice of the lattice action and the lattice definition of the topological 
charge. Although on the lattice the action barriers are typically finite, they may 
cause problems in the numerical sampling of the various Q sectors. In particular, a 
standard algorithm (e.g. the hybrid Monte Carlo algorithm) that changes the gauge 
field configuration in small steps may get stuck in a fixed topological charge sector 
because it cannot overcome the action barriers. The question arises how one can 
extract meaningful physical information from numerical calculations that get stuck 
in a given topological sector. Locality and cluster decomposition properties sug- 
gest that the restriction to fixed topology becomes irrelevant in the infinite volume 
limit. Still, in order to extrapolate reliably to that limit, one must understand the 
corresponding finite size effects. 

The Atiyah-Singer index theorem, Q — ul — nu, connects the physics of light 
quarks to the topology of the gluon field: the difference of the number of left- and 
right-handed zero modes of the Dirac operator, Ul and nji, is equal to the topological 
charge. Consequently, when one quark is exactly massless, the zero modes of the 
Dirac operator eliminate all nontrivial topological charge sectors (with Q ^ 0) from 
the QCD partition function. Hence, working at fixed Q = is then correct even 
in a finite volume. However, even in that case various physical observables (e.g. 
the chiral condensate) get contributions from Q ^ sectors, and one still needs 
to sample several topological sectors in order to compute those quantities. In the 
presence of a massless quark, the complex phase of the Boltzmann factor exp{i9Q) 
for a ^-vacuum is trivial and the physics is independent of 9. When all quarks are 
massive, there are nontrivial ^- vacua effects and all topological sectors contribute to 
the partition function. In small space-time volumes with f3V{^'^)m <^ 1 (where f3 
is the inverse temperature, V is the spatial volume, (^f^') is the chiral condensate, 
and m is the quark mass) the partition function is still dominated by the Q = 
contribution. On the other hand, in large volumes with pV{^^)m ^ 1 one must 
sample a large number of topological charge sectors. 
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When standard lattice fermion actions are used, the relation between the topol- 
ogy of the gluon field and the zero modes of the Dirac operator is obscured by lattice 
artifacts. In that case, the effects that we discuss here are difficult to observe in a 
numerical computation. However, when one uses a fermion action that obeys the 
Ginsparg- Wilson relation Ij, the chiral properties of continuum fermions persist on 
the lattice. Our results — which we derive for the continuum theory — directly 
apply to this type of lattice fermions. 

The rest of the paper is organized as follows. In section 2 the basic physics of 
^-vacua and the topological charge is reviewed, and in section 3 a simple formula for 
the fixed topology finite size effect on particle masses is derived. In section 4 chiral 
perturbation theory is used to study the ^-dependence of the vacuum energy and the 
pion mass. In section 5, we explore the large Nc limit in chiral perturbation theory. 
Since current lattice calculations are performed far from the chiral limit, in section 
6 the Q-dependence of the T^'-mass arising in an instanton gas model is estimated 
and compared with lattice results. Finally, section 7 contains the conclusions. 



2 Topological Charge Sectors and 6^- Vacua 

The partition function of QCD in a ^-vacuum is given by 

Z{e) = J VAV^V^ exp{-S[A,^,^])exp{-ieQ[A]), (2.1) 

where 

= ^/ d'x e^,p,TrF^,Fp, (2.2) 

is the topological charge of the gluon field A^(a;) = zA^(a;)A" with field strength 

F^,{x) = df,A,{x) - d,Af,{x) + [A^{x), A,{x)], (2.3) 

and S[A, ^, \E'] is the Euclidean QCD action, which is invariant under gauge trans- 
formations 

A'^ix) = g{x){A^{x) + d,)g{xy, ^'{x) = g{x)^{x), ^'{x) = ^{x)g{xy. (2.4) 

We consider the theory in a finite spatial volume V = L^LyL^ with periodic bound- 
ary conditions at temperature T, i.e. with finite Euclidean time extent Lt = (3 = 
1/T. In order to allow for nontrivial topology, it is necessary to impose periodicity 
only for gauge invariant quantities. Thus, A^{x), \E'(x) and ^{x) are periodic only 
up to gauge transformations 

A^{x + L^fx) = gf,{x){A^{x) + d^)g^{x)\ 

^{x + L^fi) = g^{x)'^{x), ^(x + L^/i) = ^{x)g^,{x)\ (2.5) 
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Here fi is the unit- vector in the /^-direction and the gfj.{x) are a set of transition 
functions that define the boundary condition. As first reahzed by 't Hooft j2], for 
an SU{Nc) Yang-Mills theory these boundary conditions fall into gauge equivalence 
classes labeled by a twist tensor that takes values in the center of the gauge group. 
In general, for a nontrivial twist tensor, the topological charge can take fractional 
values that are multiples of 1 /iVc In the presence of matter fields that trans- 
form nontrivially under the center (e.g. quarks in the fundamental representation of 
SU{3)), single valuedness of the matter field implies that the twist tensor is trivial. 
Hence, in the QCD case the topological charge is an integer and not a fraction. It 
is interesting to note that in supersymmetric Yang-Mills theory the gluinos, which 
transform under the adjoint representation, leave the center symmetry unbroken. 
Consequently, one then expects fractional topological charges. The role of the topo- 
logical charge in QCD at finite volume has been worked out in great detail by 
Leutwyler and Smilga |3]. Although eq. ()2.ip is quite formal, it can be given a con- 
crete meaning on a space-time lattice using, for example, Ginsparg- Wilson fermions. 
Below we will ignore this subtlety and perform formal manipulations on eq. (j2.H) . 
assuming that all of these can be given a concrete realization on the lattice. 

If we define H to be the QCD Hamiltonian with the energy eigenstates \n, 6) in 
a given ^-vacuum, then we can write 

Z{d) = Y.{n,e\eM~l3H)\n,e) = Y.ei<.v{-f3Er,{e))- (2-6) 

n n 

For sufficiently large spatial volumes, the 6*- vacuum energy is given by Eq{6) = 
Veo{6), where 60(6*) is the vacuum energy density. For small values of 6 it is given 
by 

eo{e) = ^XtO' + lO' + -, (2.7) 

where 

Xt - (2.8) 

is the topological susceptibility (at 9 = 0). For P[Ei{6) — Eo{6)] ^ 1 the vacuum 
state dominates and we get 

Zi9) = expi-pVeom = exp{-^9')[l + 0{j9% (2.9) 

For small values of 6 and for M^r-L ^ 1 (where L = = Ly = L^) the first excited 
state is the pion, i.e. Ei — Eq = M^^. On the other hand, in the limit of massless 
quarks chiral symmetry is restored in a finite volume and a rotor spectrum of states 
exists above the vacuum In that case, the energy of the first excited state is 
given by 

^i(0)-^o(0) = ^^, (2.10) 

where Nj is the number of quark fiavors and is the pion decay constant. In 
other words, the energy gap is inversely proportional to the spatial volume. As the 
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quark mass increases, this state evolves into the state of a pion at rest. For 6 = n 
(depending on the quark masses and Nf) the CP symmetry breaks spontaneously 
In that case, there are two degenerate 6'- vacua in the infinite volume. In a finite 
volume, due to tunneling, there is an exponentially small energy gap 



Ei(7r) - Eo{7i) =Aexp{-aV). 



(2.11) 



Here a is the tension of a 3-dimensional interface separating Euclidean space-time 
regions that are in distinct CP phases. 



The QCD partition function at fixed topological charge Q is given by 
Zq = J VAV^V^ 6Q,Q[A]exp{-S[A,^,^). 
Writing the Kronecker (5-function as 

^Q,QIA] = ^ exp[teiQ - Q[A])], 

the fixed Q partition function turns into an integral over all 6 values 
Zq = — d9 Z{9) expitOQ) 



(2.12) 



(2.13) 



1 

2'n J-n 
1 



27V 



(2.14) 



At sufficiently low temperatures (/5i?i(0) ^ 1) and sufficiently small topological 
charges (i.e. \Q\/\/]3V remains finite) the partition function reduces to 



1 

2^ 



de exp{-pVeo{9))exp{i9Q). 



(2.15) 



This equation is the starting point to derive the effects that we are interested in. 
For large space-time volumes jSV, the fixed Q partition function is dominated by 
the ^-vacua with a small energy density. Since the 6 = vacuum has the smallest 
energy, we see that the physics at large PV is dominated by 6* = 0. Hence, we can 
use a saddle point approximation to obtain a formula for Zq. We assume a large 
PV keeping Qj^f^V finite and fixed. Using eq. ()2.7j) one obtains 



I 2tx 

W) 



exp 



exp 



2(3Vxt 

2m 



i + o( 
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/3V' 



(2.16) 
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3 Effective Mass Determination at fixed Q 



It is possible to determine the hadron spectrum from a Monte Carlo simulation even 
if it gets stuck in a fixed topological sector. Typically hadron masses are extracted 
from 2-point correlation functions of operators with the appropriate quantum num- 
bers. Higher n-point functions are important in the determination of other physical 
quantities like matrix elements and coupling constants. In a 6'-vacuum a general 
n-point correlation function of operators Oi takes the form 

= Y{f)j Or02...0neM-S[A,^M)^M-^SQ[A]). (3.1) 

Hence, in a sector of fixed topological charge one obtains 

1 1 



Gq = {o,02...0r:)Q = r de z{e)G{e) exp(^eg). (3.2) 



Again, assuming small enough temperatures (/5i?i(0) ^ 1) and sufficiently small 
topological charges (i.e. a finite \Q\/ ^/JW)^ Zq is given by eq. ()2.16|l and Z{0) is given 
by eq. ()2.9j) . Assuming that {0i02---0n)e is smooth over the range 6*^ < 1/ [3Vxt = 
1/{Q'^), one can again use the saddle point method to evaluate the integral over 6. 
Then one obtains 

GQ = Gie,) + ^^G"ies) + ..., (3.3) 

where a prime denotes a derivative with respect to 6. The (purely imaginary) value 
of 9 at the saddle point is given by 

Let us now consider a 2-point correlation function 

^ ''VAV^V^ 0{ti)0{t2)exp{-S[A,^,^])exp{-teQ[A]) 



z{e) 

1 

m 



m 

\- Y,{n, e\Oexp{-H{h - t2))0 exp{-H{f3 - ti + ts)) |n, 6) 

^Y.^M-f3E^{omm,e\o\n,e)\^exp[-{En{e)-E^{e)){h-t2)], 

(3.5) 



relevant for determining particle masses. Here the operator O has the appropriate 
quantum numbers in order to create the particle of interest at rest from the 6 = 
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vacuum. The 6'-dependent particle mass is then given by M{6) = Ei{6) — Eq{6). 
For large time separations (ti — and at sufficiently low temperatures (large /3) 
the correlation function reduces to 



{Oih)Oih))9 = \{0,e\O\l,e)\^exp[-Mi9)ih~h)]. 



(3.6) 



Using eq. ()3.3|l and assuming 




(3.7) 



one obtains 



(3.8) 



where Aq is a constant and 




) + ... 



f^Vxt 




) + .... 



(3.9) 



It should be noted that the concept of an effective "mass" Mq (which is independent 
of ti — ^2)) in a fixed topological sector makes sense only to the quadratic order of the 
saddle point expansion considered here. Higher order terms lead to time- dependent 
"masses". This is not surprising because, due to the topological constraint, QCD at 
fixed Q is a nonlocal theory that does not even have a Hamiltonian. We note again 
that we take the large PV limit keeping Qj ^ fiV fixed. Hence, the term Q^/ (Q^) is 
of order 1. When Q itself is fixed this term becomes of order Xj^V . In this case, 
there are other terms of the same order that we have not calculated. Hence, in order 
to be consistent, this term should then be dropped. 

The above discussion shows that if a calculation performed in a space-time vol- 
ume PV gets stuck in a fixed Q sector, one would extract a mass Mq that deviates 
from the infinite volume result by a term proportional to l/(/9V"). This term de- 
creases rapidly as one increases the space-time volume. Of course, from a simulation 
that is stuck in a fixed topological charge sector, one will not be able to extract (Q^), 
because it will be impossible to determine the relative weight of different sectors. 
Still, one should be able to perform separate computations for the various topologi- 
cal sectors, simply by starting the calculation in a given Q sector. Then, in principle, 
one can measure Mq in the various sectors and at various space-time volumes and 
thereby extract M(0) , M"(0) and (Q^) using eq. (|3.9p . It is interesting to note that 
when Mq is averaged over Q with the distribution exp(— (5^/2((5^)) one obtains the 
mass M(0) in the 6' = vacuum. 

Let us briefly consider the extreme case of zero quark mass. Then the pion is 
massless and the box size is always small compared to the infinite Goldstone boson 
correlation length. For m = there are no ^-effects — ^(^) is independent of Q 
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— and the fixed Q partition function Zq is non-zero only for (5 = 0. Hence, in 
tliis case, it is correct to work at fixed topological charge (5 = 0, unless one wants 
to measure observables like the chiral condensate that receive contributions from 
(5 7^ sectors. 



4 ^-Dependence of Physical Quantities 

Based on the previous discussion, it is clear that in a large space-time volume the 
effects of fixed topology can be determined from the 6'-dependence of certain physical 
quantities. For example, the partition function Zq at fixed topological charge Q can 
be determined from the 6'-dependence of the vacuum energy 60(6*) using eq. ()2.15|) 
and the effective mass of a hadron Mq can be obtained from M{9) using eq. ()3.9p . 

Until now we have written the QCD action in a 6'- vacuum as S[A, \E', ^] +i6Q[A\. 
It is well known that, using an anomalous U{1)a chiral transformation, the 6- 
dependence can be moved into the quark mass matrix M. = diag (m^, rrifi, ...,miy^) 
(for Nf flavors), i.e. 

9 •' 

+ j (fx {^LMe'''''f<^R + ^RM''e~'''''^^L}. (4.1) 

Thus, 9 appears simply as the complex phase of the determinant of the mass matrix 

Although 9 enters the QCD action in a simple way, in general its effects on 
physical quantities are difficult to determine. However, for small quark masses chiral 
perturbation theory allows us to understand such effects. It is well known that the 
low-energy physics of Goldstone bosons in QCD can be described systematically 
using an effective theory. In the case of Nf light quarks to lowest order the action 
of the Goldstone boson field U e SU (Nf) takes the form 

d^x {^Ti[d^U^d^U] - A^Tr[A^e'^/^^t/t + M^-^'^^'fU]}, (4.2) 
4 ^N f 

where F^^ is the pion decay constant. The term proportional to F^ is invariant under 
chiral SU{Nf)L ® SU{Nf)R transformations 

U'{x) = LU{x)R\ (4.3) 

(with L G SU{Nf)L and R G SU{Nf)ji) while the mass term explicitly breaks chiral 
symmetry. 
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First, we consider the case of Nf degenerate flavors of mass m. Then the mass 
matrix takes the form Ad = ml. In this case, the chiral symmetry is exphcitly 
broken down to the vector subgroup SU{Nf)L=ji. When the vacuum angle changes 
by 27rn, exp{i6/Nf) changes by an element z = exp{27ini/Nf) G ^.{Nf) — the center 
of the flavor group SU{Nf). Although the action itself is not 27r-periodic in 6, the 
corresponding partition function 

Z{e) = jvU exp(-5[f/]) (4.4) 

is, because the shift exp{i9' /Np) = z exp{i9/Nf) can be compensated by a redefi- 
nition of the field U' = zU. For 6 G [— 7r,7r], the effective action is minimized by 
constant field configurations f/ = 1. The value of the action at the minimum 

S[l] = -(3V{^^)mcos{e/Nf) (4.5) 

determines the ^-vacuum energy density (normalized to zero at 6' = 0) as 

eo{e) = {^^)m{l - cos{e/Nf)). (4.6) 

This expression is valid for 9 G [—27i/Nf,27i/Nf]. When extended periodically to 
other ^-values, it has cusps at ^ = 47r(n+ \)/Nf indicating a phase transition in the 
vacuum angle 0. Comparing this result with eq. ()2.7|) . we identify the topological 
susceptibility as 

X. = (4.7) 

In order to determine the ^-dependence of the Goldstone boson mass, we expand 
U = exp^iii"' X"' / Fj,) in powers of vr" and read off M^{9) as the coefficient of the tt^tt'^ 
term. This yields 

M2(^) = M^2(0)cos(^/iV^), (4.8) 

where 

Mm = -j^, (4.9) 

is the Goldstone boson mass squared at 6* = 0. 

Using leading order baryon chiral perturbation theory [7] the quark mass depen- 
dence of the nucleon mass (at 6' = 0) is given by 

M40) = M°[l + c^]-|i^, (4.10) 

where and qa are the nucleon mass and the neutron decay constant in the chiral 
limit, and c is a known combination of low-energy constants. At this order in chiral 
perturbation theory the ^-dependence enters only through the pion mass of eq. (|4.8p 
such that 
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Let us now consider the case of two non-degenerate flavors with masses mu 
and m^. Then the mass matrix takes the form Ai = diag(m„,md), and chiral 
symmetry is exphcitly broken down to U{1)l=r. Now the action is minimized for 
U = dia.g{exp {iip),exp{—iip)) where 



tamp 



rriu - rrid 
ITLu + nid 



tan ■ 



e 



(4.12) 



and one obtains 



rriu + rud 



1 — cos 



(m„ - rUdY , 2 ^ 

1 + 7 TTT tan"^ - 

[niu + rridY 2 



(4.13) 



The pion mass is given by 



Ml{d) = Ml{Q)co^ 



2 y ^ (m„ + rridy 



tan 



(4.14) 



where now 



2,^, (^^)(m„ + md) 



2F2 



(4.15) 



5 ^-Dependence at large Nc 

The flavor-singlet pseudoscalar particle 1]' is special in QCD. Due to the anomaly it 
is heavy (and thus not a Goldstone boson) and hence its physics cannot be under- 
stood using chiral perturbation theory. However, in the large Nc limit the anomaly 
disappears and the t]' becomes a Goldstone boson which can now be studied in chiral 
perturbation theory. The low-energy effective action of eq. fl4.2|l is modified to jBj 

4 ^N f 

+ ^x?(^logdetf/)2. (5.1) 



Here [/ is a U{Nf) matrix with complex determinant exp(i^2Nfri' / F^^). The last 
term in eq. ()5.1|) gives rise to the rj' mass 

M,^, = (5.2) 

At large N^, the pion decay constant goes as ~ Nc while the topological suscep- 
tibility Xt of the quenched (zero flavor) theory is of order 1. Hence, as Nc goes to 
infinity, the r/'-meson becomes a massless Goldstone boson. The occurrence of the 
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rl field in ttie chiral Lagrangian tias a dramatic effect on the ^-dependence. Since Q 
enters as the determinant of the mass matrix, it can be absorbed into a redefinition 
of the field Vj such that 

S\IJ\ = [d^x{^Tr[d^U^d,,U]-^Ti[MU^ + M^U]} 

+ ^Xti^^ogdetU-9f. (5.3) 

At large Nc the last term in eq. ()5.3p can be neglected and thus all ^-dependence 
disappears from the theory even at non-zero quark masses. This is consistent because 
the anomaly disappears in the large Nc limit and Ua{^) chiral rotations can be used 
to relate different 6'-vacua which thus become physically indistinguishable. 

At large but finite Nc we can use the action of eq. ()5.3|) to derive the 6'-dependence 
of the Goldstone boson masses for Nf degenerate quark flavors of mass m. First, 
we must find the minimum of the action, which turns out to be at 



U = exp{z^2/Nfr]',{d)/F^)l, (5.4) 

where 

mmM^^^-^) + x'tNfi^f"^ -e) = 0. (5.5) 

When one is closer to the chiral than to the large Nc limit the second term on the 
left-hand side of eq. fl5.5p dominates and then 



It is easy to show that the vacuum energy takes the form 

eo(^) = -mm) cos(^^) + ^X?(V^^ - Of. (5.7) 
The ^-dependence of the pion mass 

Ml{e) = Ml{Q)cos{^^-^\ (5.8) 



is similar to the one in eq. ()4.8|) except that 6 is now replaced by ^j2NfriQ{9)/ F^. 
The 7]' mass is given by 

Mi{e) = ^-^ + Mm- (5.9) 

IT 

At Nc = oo the first term on the right-hand side of eq. (j5.9p vanishes, thus making 
rj' degenerate with the pion. Solving explicitly for ?7q, we also note that in the limits 
m ^ or A'^c ^ oo physical quantities become 6'-independent. 
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Recently, the dependence of M^r and M^/ on Q has been studied in lattice calcu- 
lations [S]. Each mass was evaluated in two domains, one with IQI < 1 and the other 
with \Q\ > 2. While the pion masses in both domains agreed with each other within 
statistical errors of a few percent, the r/'-mass was about 15 percent heavier in the 
|Q| > 2 domain than in the IQI < 1 domain. Unfortunately, the quark masses used 
in current lattice calculations are far from the chiral limit such that the above results 
from chiral perturbation theory should not be applicable to these lattice data. For 
example, in a typical full QCD lattice calculation Mj^/Mp > 0.6. Still, it is tempting 
to compare our predictions with the lattice data. Assuming that we are closer to 
the chiral than to the large Nc limit, using eqs. (j5.5l3.9|) we find 



M^Q = M,(0) 



1-^^(1-^) 



(5.10) 



In the lattice calculation of ,8j I3V = 6.81fm'^, Xt = 0.70fm ^, and hence (Q^) = 
PVxt = 4.75. This means that the predicted shift of the effective pion mass at Q = 
is given by -M^{0)/{ANf{Q^)) = 0.013M^(0). Hence, the effect of fixed topology 
is of the order of 1 percent, consistent with the lattice results. 

At large but finite A^^^ and close enough to the chiral limit (i.e. when eqs. ()5.8l5.9|) 
apply) the second 6'-derivatives M^'(O) and M"{0) are equal. Hence, using eq. (j3.9|) 
we would expect that the Q-dependent shift of the effective masses of rj' and the 
pion are the same, in disagreement with the lattice results. 



6 Q-Dependence in an Instanton Gas 

It is useful to try to understand how fixed topology affects QCD in the domain ex- 
plored by actual lattice calculations. Motivated by the phenomenological success of 
the Veneziano-Witten formula relating the ?7'-mass to the quenched topolog- 
ical susceptibility Xt , in this section we will seek a qualitative understanding of the 
Q-dependence of the r^'-mass by modeling the fluctuations of the topological charge 
by an instanton gas. 

We begin by assuming independent Poisson distributions for instantons and anti- 
instantons in a space-time volume (3V with (A^) = (A") = A. Then the probability 
of having A^ instantons and A^ anti-instantons is given by 

\iN+N) 

Using Q = N — N one obtains 

(^Q^) = (^N + N) = 2X = f3Vx't. (6.2) 



12 



Fixing the topological charge to Q the probability distribution takes the form 

1 1 r A^+^ 



where 



^Q(A^,Ar) = — -y_ de j^e-'^ - {N - N))l (6.3) 

Zq = ^ r de exp[-2A(l - cos^)] exp(i^Q). (6.4) 



Motivated by eq. ()6.2p we want to model the Q-dependence of the T^'-mass by 
assuming 

_2Nf {N + N)Q 
In the large 13V limit (with Q/ \/(3V fixed) one obtains 



(iV + iV)g_^0 1 



1 - 



2 



(6.6) 



This is a simple and physically appealing result. When M'^,q is averaged with 
the distribution Zq ~ exp{—Q^/2{Q'^)) in the appropriate limit, the r]' mass is 
consistent with the Veneziano-Witten formula. The finite volume corrections provide 
the desired Q-dependence. Using eqs. fl6.5l6.6p for M^, , we find that the percentage 
shift of M^/ at Q = is given by l/4(3Vxt = V4(Q^)- In the lattice calculation 
of [S] 1/4(Q2) ^ 0.06, which is of the same order of magnitude as the observed 10 
percent shift in the 77'- mass between the sum over all Q and the sum over \Q\ < 1.5. 



7 Conclusions 



In this work, we have shown how to correct finite size effects that occur in lattice 
QCD calculations at fixed topology. We have derived a formula that relates the 
Q-dependent effective mass Mq to the true mass M(0) and its second derivative 
M"(0) at vacuum angle ^ = 0. The difference Mq — M(0) is of the order 1/pV and 
thus vanishes quickly in the large volume limit. By observing the predicted finite 
size effect in a lattice calculation that gets stuck in fixed topological charge sector, 
one can identify not only M(0) but also M"(0) and thus one can learn something 
about the 6'-dependence. Close to the chiral limit, we have used chiral perturbation 
theory to predict M"(0) for the pion and the nucleon. For large we have also 
derived this quantity for the r^'-meson. Our formulae are in reasonable agreement 
with existing lattice data for the effective pion mass at fixed Q, despite the fact 
that those calculations are performed far from the chiral limit. On the other hand, 
perhaps not surprisingly, the observed Q-dependent effective T^'-mass does not agree 
with the large chiral prediction. Interestingly, a dilute instanton gas model is 
qualitatively consistent with the lattice data. 
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